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Abstract 

Three kinds of the response mechanisms to the external pressure have been found for double- 
walled carbon nanotube (DWNT) bundle, depending strongly on their average radius and sym- 
metry. The small-diameter DWNT bundle undergoes a small discontinuous volume change, and 
then deform continuously. The intermediate-diameter DWNT bundle collapses completely after 
a structure phase transition (SPT). Significantly, two SPTs exist for the larger-diameter DWNT 
bundle if the outer tube has no Cq or C3 symmetry. It would be interesting to search for signatures 
of these different structural transformations by experimentally investigating mechanical, optical 
and thermal response functions of DWNT bundle. 
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It is well known that the physical properties of the carbon nanotubes (CNTs) depend 
much on their geometrical structures, and so can be easily changed by an applied pressure 
or strain, which could be used to fabricate the nanoscale electromechanical coupling devices 
and transducers. For example, a uniaxial strain on the single-walled carbon nanotubes 
(SWNTs) can cause a metal-semiconductor transition^. 

Recently, the effect of hydrostatic pressure on the CNT bundle has attracted much 
attention in experiments^'^'^'^'^i^'^i^'i^, including the Raman spectroscopy, and theoretical 
calculations^'i^ii^'^. The study of the SWNT bundles indicates the Raman peaks shift 
to higher frequencies with increasing hydrostatic pressure, and the radial breathing mode 
(RBM) disappears from the spectrum above the critical pressure, showing a SPT. On the 
other hand, there are much fewer studies of the hydrostatic pressure's effect on the DWNT 
bundle, which is the simplest multi- walled carbon nanotubes (MWNTs). Very recently, the 
Raman measurements on the DWNT materials under hydrostatic pressure^i^i^ have found 
the RBM intensity of the outer tubes decreases rapidly with increasing pressure, exhibiting 
a similar behavior to that of SWNT-i^i^i^i^iiiSiiiiii^ii^, but the inner tubes appear to be con- 
siderably less sensitive to the pressure, e.g., the pressure coefficient of the inner tube is 45% 
smaller than that of the outer tube. 

Therefore, in this work, we intend to study the structural change of DWNT bundles 
under hydrostatic pressure. Our numerical results show that the symmetry and diameter 
of outer tube decide mainly the change of cross section and response behavior of DWNT 
bundles under the pressure. More importantly, two SPTs are found to exist in some DWNT 
bundles if their outer tubes have the larger diameters and no Cq or C3 symmetry. 

The zero-temperature structural minimizations of the enthalpy {H = U + PV) were 
carried out on a supercell containing a 2 x 2 x 2 DWNT bundle using the universal force 
field (UFF) methodr^i^. In order to induce the SPT, a step-wise increasing hydrostatic 
pressure was applied to the DWNT bundle, minimizing the enthalpy of the DWNT bundle 
after each pressure increment. 

For (5,5)@(10,10), (7,7)@(12,12) and (9,0)@(18,0) DWNT bundles, we also perform the 
first principles structure relaxations using the total energy plane-wave potential method - 
in the framework of local density approximation (LDA). The ion-electron interaction was 
modeled by the projector augmented wave (PAW) method^ with a uniform energy cutoff of 
500 eV. The smearing width was taken to be 0.04 eV in the ground state. And a very good 
agreement is obtained between the force field and the first principles calculation results (see 
Table I), indicating that the UFF method is suitable for the system researched in this paper. 

The transition pressures, Pc and P^, for different DWNT bundles are listed in Table II. 
As we shall see, there exist three kinds of the response behaviors to external pressure, and 
their loading curves vs hydrostatic pressure are given in Fig. 1(a)- 1(c), respectively. 

Firstly, for small-diameter (5,5)@(10,10) DWNT bundle, a small discontinuous volume 
change appears at P = 18.01 GPa, accompanied by a cross section's change between two 
deformed hexagons, as seen from Fig. 1(a). At the same time, the inter tube (5,5) could 
not collapse with a minimum distance of 4.35 A between its two opposite walls, which is 
larger than the distance of 3.4 A between nearby layers in the turbostratic graphite. If the 
pressure is increased further, this distance will approach continuously to 3.4 A. On the other 
hand, our simulation indicates that the loading curve of (5,5) SWNT bundle [also shown 
in Fig. 1(a)] varied with external pressure is continuous too and no obvious SPT happens. 
But (10,10) SWNT bundle collapses at P^ = 3 GPa, forming a peanut-shaped cross section 
with a separation of about 3.4 A between its two opposite parallel walls. So, it is clear that 
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existence of an inner (5,5) tube increases the ability of the outer (10,10) tube to resist the 
apphed pressure, making so the (5,5)@(10,10) DWNT bundle do not collapse at 18.01 GPa. 

As for the intermediate-diameter DWNT bundles, e.g., the (7,7)@(12,12), (9,9)@(14,14), 
(10,10)@(15,15), (11,11)0(16,16), (12,12)0(17,17) and (13,13)0(18,18), as shown in Fig. 
1(b), it is found that all of them undergo one SPT and collapse completely. Our simulations 
reveal similar collapses exist for the (7,7), (9,9), (10,10), (11,11), (12,12) and (13,13) SWNT 
bundles. Taking (10,10)0(15,15) DWNT bundle as an example, its collapse pressure Pd = 
4.68 GPa that is higher than either that of (10,10) SWNT bundle {Pd = 3 GPa) or that of 
(15,15) SWNT bundle {Pd = 1.3 GPa). This means that the outer tube acts as a protection 
shield, and the inner tube supports the outer one and increase its structure stability. Our 
results are consistent with the experimental measure results^i'^'^. 

However, the response behaviors of the larger-diameter DWNT bundles become complex, 
as shown in Fig. 1(c). The (16,16)0(21,21) and (19,19)0(24,24) DWNT bundles still 
collapse after one SPT. Surprisingly, the (15,15)0(20,20), (17,17)0(22,22), (18,18)0(23,23) 
and (21,21)0(26,26) DWNT bundles undergo two different types of SPT with the transition 
pressure Pc and Pd respectively, accompanied by their cross sections change from deformed 
hexagon to racetrack first and then to peanut shape, and the DWNT bundles collapse 
completely after the second SPT. It should be noted that the DWNT bundles exhibit some 
small discontinuous volume change between the different racetrack-shaped cross sections. 
We think two SPTs and the racetrack-shaped cross sections could be observed in future 
experiments if the pressure increment is selected properly. 

Furthermore, it is interesting to ask why the racetrack-shaped cross sections do not 
appear for the (16,16)0(21,21) and (19,19)0(24,24) DWNT bundles before their final col- 
lapse? Our simulations reveal that not only their cross sections, but also the ones of the 
intermediate-diameter DWNT bundles: (7,7)0(12,12), (10,10)0(15,15), (13,13)0(18,18) re- 
main relatively the better hexagon with almost equal sides and corners just before the SPTs, 
which are different from the other DWNT bundles with the deformed hexagonal or race- 
track cross sections before the SPT (in Fig. 1(b)- 1(c) we only show the cross sections of 
some DWNTs, e.g., the (9,9)0(14,14), (10,10)0(15,15), (13,13)0(18,18), (16,16)0(21,21) 
and (18,18)0(23,23) as examples). It is known that the symmetry group of the isolate (n,n) 
and (n,0) SWNT is T^Dnh- The symmetries of and ah are retained under pressure, but 
in the bundle the cross section symmetry is reduced from Dn to Cg for the specific outer 
tubes: (12,12), (18,18) and (24,24), and to C3 for the (15,15) and (21,21) tubes. A neces- 
sary condition to form an ideal hexagonal lattice is that the tube itself has a rotational 
axis, and C3 symmetry can also be somewhat matched to the hexagonal lattice symmetry. 
Thus, with increasing pressure, these outer tube's cross section can keep better hexagon 
before the SPT. More importantly these DWNT bundles only undergo one SPT to reach 
the stable collapsed structure no matter how big is the diameter. On the other hand, the 
DWNT bundles in which only their inner tubes have the Cg or C3 symmetry, for example 
the (9,9)0(14,14), (12,12)0(17,17), (15,15)0(20,20), (18,18)0(23,23) and (21,21)0(26,26), 
have still the deformed hexagonal or racetrack-shaped cross sections just before the SPT, 
and two SPTs would happen with increased diameters. So, the symmetry of the outer tube 
is a very important factor to decide the cross section shape of the DWNT under the pressure, 
and comparatively the effect of their inner tube is very small. 

We have also studied the structure transformations of zigzagOzigzag DWNT bundles 
under the hydrostatic pressure, which are found to show a qualitatively same response 
mechanism to external pressure with the armchairOarmchair DWNT bundles, as indicated 
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by Table 11. 

Finally, in order to have an insight into the relationship of the collapse pressure Pd with 
the DWNT diameter and its tube symmetry, we introduce the DWNT's average radius 
Rave, which is defined as an average value of the outer and inner tube's radius of DWNT. 
Here we show the variation of Pd with Rave in Fig. 2, from which it can be seen that the 
collapse pressure can be well fitted to ~ l/-RLe DWNT bundles. It is known that 

the tube-tube coupling in a bundle is described by the van der Waals force, which has a 
small effect on the bundle's collapse pressure. And both the SWNT and DWNT can be 
described by the continuum elasticity theory as the continuous hollow cylinders. So their 
bundles would have similar response to the hydrostatic pressure, and the only difference is 
that the effective radius of DWNT is the average of inner and outer tubes' radius. On the 
other hand, it has been proved by the continuum elasticity theory [19] and the molecular 
dynamics simulations [20] that the Pd of an individual SWNT is inversely proportional 
to its cubic radius. Thus, it is reasonable that the relation of Pd ~ l/-RLe appears in 
the DWNT bundle. In addition, for the DWNT bundles in which the outer tube has Cq 
or C3 symmetry, the collapse pressure is larger than the fitted result. This phenomenon 
is more obvious, especially for the zigzag@zigzag DWNT bundles, e.g., the (9,0)@(18,0), 
(15,0)@(24,0), (21,0)@(30,0), (27,0)@(36,0) and (33,0)@(42,0). The reason is that both of 
their outer and inner tubes have Cq or C3 symmetry, further enhancing the matching to 
the hexagonal lattice and increasing the structure stability of the system. This means the 
matching between the DWNT symmetry and the lattice symmetry can increase the ability 
of the DWNT bundle to resist the applied pressure. 

In addition, we have also made simulation on the (11,2)@(12,12) DWNT bundle, which 
is composed of two coaxal SWNTs with different chiral angles, and the result is shown in 
Table I, Fig. 1(b) and 2. The period of tube (11,2) is just seven times larger than that 
of tube (7,7), and its radius is equal to that of (7,7). It can be found from Fig. 1(b) that 
the (7,7)@(12,12) and (11,2)0(12,12) DWNT bundles have the slightly different collapsed 
pressures of Pd = 10.6 GPa and 11.51 GPa, respectively, i.e. the chiral symmetry of the 
inner tubes in the DWNT bundle has a smaller effect on their collapse pressures. 

In summary, our calculations show that the structural transformations of DWNT bun- 
dles under hydrostatic pressure is different from those of the SWNT bundles^'^'^'^i^'i^i'ii'iSi'i^ 
and isolate DWNT—. One or two SPTs exist depending on the symmetry and diameter 
of DWNT bundles, which manifest the complexity of nanotubes. It would be interesting 
to experimentally determine mechanical (e.g., compressibility), optical (e.g., Raman spec- 
trum) and thermal (e.g., heat capacity) response functions of DWNT bundles to search for 
signatures of these different types of structural transformations. 
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TABLE 



TABLE I: Calculated lattice parameters of some DWNT bundle without external pressure by UFF 
method and first principles method, a, b and c are the lattice constants of DWNT bundles, and 
a, 13 and 7 are the angles between the two lattice vectors. 
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TABLE II: Calculated critical transition pressure Pc and Pd of DWNT bundles. Rave is the average 
value of the inner and outer tube's radius in the isolate DWNT directly rolled up from the ideal 
graphene sheet. The suffix (NC) means at this pressure, no collapse happens. 
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FIG. 1: (Color online) Loading curves for different armchair@armchair DWNT bundles as a 
function of hydrostatic pressure, (a) small-diameter (5,5)@(10,10) DWNT bundle; (b) some 
intermediate-diameter DWNT bundles; (c) the larger-diameter DWNT bundles. 
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FIG. 2: (Color online) Collapse pressure Pd as a function of the average radius of DWNT. can 
be well fitted to ~ iZ-^Le- 
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